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ABSTRACT 

We investigate the afterglow properties and large-scale environments of several short-duration gamma-ray 
bursts (GRBs) with sub-arcsecond optical afterglow positions but no bright coincident host galaxies. The pur- 
pose of this joint study is to robustly assess the possibility of significant offsets, a hallmark of the compact object 
binary merger model. Five such events exist in the current sample of 20 short bursts with optical afterglows, 
and we find that their optical, X-ray, and 7-ray emission are systematically fainter. These differences may be 
due to lower circumburst densities (by about an order of magnitude), to higher redshifts (by Az « 0.5 - 1), or to 
lower energies (by about a factor of 3), although in the standard GRB model the smaller 7-ray fiuences cannot 
be explained by lower densities. To study the large-scale environments we use deep optical observations to 
place limits on underlying hosts and to determine probabilities of chance coincidence for galaxies near each 
burst. In 4 of the 5 cases the lowest probabilities of chance coincidence (P(< SR) ~ 0.1) are associated with 
bright galaxies at separations of SR ~ 10", while somewhat higher probabilities of chance coincidence are asso- 
ciated with faint galaxies at separations of ~ 2". By measuring redshifts for the brighter galaxies in three cases 
(z = 0.111, 0.473, 0.403) we find physical offsets of « 30-75 kpc, while for the faint hosts the assumption of 
Z > 1 leads to offsets of ~ 15 kpc. Alternatively, the limits at the burst positions (> 26 mag) can be explained 
by typical short GRB host galaxies (L m 0.1 - 1 L*) at z > 2. Thus, two possibilities exist: (i) ~ 1/4 of short 
GRBs explode ~ 50 kpc or ~ 15 kpc from the centers of z ~ 0.3 or z > 1 galaxies, respectively, and have fainter 
afterglows due to the resulting lower densities; or (ii) ~ 1/4 of short GRBs occur at z > 2 and have fainter af- 
terglows due to their higher redshifts. The high redshift scenario leads to a bimodal redshift distribution, with 
peaks at z ~ 0.5 and z ~ 2. The large offset scenario leads to an offset distribution that is well-matched by 
theoretical predictions of NS-NS/NS-BH binary kicks, or by a hybrid population with globular cluster NS-NS 
binaries at large offsets and primordial binaries at offsets of < 10 kpc (indicative of negligible kicks). Deeper 
constraints on any coincident galaxies to > 28 mag (using the Hubble Space Telescope) will allow us to better 
exclude the high-redshift scenario. 

Subject headings: gamma-rays:bursts 



1. INTRODUCTION 

The bimodality of gamma-ray burst (GRB) durations 
(Kouveliotou et al. 1993) is indicative of separate progenitor 
populations for long- and short-duration GRBs. While di- 
rect observational s upport exists for the mas sive star origin 
of long GRBs (e.g., Woosley & Bloom 2006), the most pop- 
ular progenitor model for short GRBs is the coalescence of 
compact object binaries with neutron star and/or black hole 
consti tuents (NS-NS /NS-BH; |FjchleretalJ U989J iPaczvnskil 
I199U iNaravan et~aT1 1 1992b . One of the key predictions of 
this model is that some systems will experience large velocity 
kicks at birth, leading to eventual mergers well outside of the 
host galaxies, in galactic halos and the intergalactic medium 
(iBloom et alJll999HFryer et al]ll999tlBelczynski et alll2006l) . 
These models predict that 10-20% of all mergers will occur 
at offsets of > 20 kpc for Milky Way mass galaxies. In such 
environments the resulting afterglow emission is expected to 
be fainter than for bursts occurring in coinc idence with their 
host g alaxies due to the low ambient density dPanaitescu et alj 
l200UlPerna & Belczvnskill2002l) . 

A subset of NS-NS binaries ( ~ 10-30%) may be f ormed 
dynamically in globular clusters (Grind lay et alJl20 06). lead- 
ing to possible large offsets and low ambient densities despite 
an absence of kicks. The distribution of offsets for such bi- 
nary systems can in principle be calculated from the spatial 

1 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cam- 
bridge, MA 02138 



distribution of globular clusters, and initial predictions are 
that 50 - 90% of these systems will h ave offsets of > 20 kpc 
(depending on the mass of the galaxy; Salva terra et al.l l2010). 
An additional expectation is that dynamically-formed bina- 
ries will be heavily skewed to lower redshifts due to the ad- 
ditional time delay bet ween the formation and core-collapse 
of the globular clusters (Hop man et alj|2005 ISalvaterra et alj 
l2008h . 

Other progenitor systems for short GRBs have also been 
proposed, including young magnetars (Thompson & Duncan 
1995), accretion-induced collapse (AIC) of neutron stars 
(Qin et al. 1998), and delayed magnetar formation through bi- 
nary white dwarf mergers or white dwarf AIC dLevan et alJ 
l2006btlMe"tzger et alj2008l) . These models are partially moti- 
vated by observations that cannot be easily accommodated in 
the standard NS-NS merger model, such as extended soft 7- 
ray emission on timescales of ~ 100 s (e.g.. | Villasenor et alj 
120051: iMetzger et all 120081; iPerlev et all l2009l) . or by phe- 
nomena such as short-duration giant flar es from soft 7- 
ray repeaters (e.g ., Hurlev et al. 2005; Pa lmer et al.1 120051 ; 
iTanvir et all 120051 : iNakar et al.l 120061) . The general expecta- 
tion is that these alternative progenitors will not experience 
kicks, and will therefore lead to bursts in coincidence with 
host galaxies. 

The detection of short GRB afterglows starting in mid-2005 
provided an opportunity to investigate the various progenitor 
models thro ugh a range of observational tests: the re dshift 
distribution dBerger et all 120071: iGal-Yam et aTJ 120081) . the 



host g alaxy de mographics (Berger 2009), t he afterglow prop- 
erties (iBergeri 12 007: Ge hrels et al.1 l2008t iKann et alJ 120081: 
iNvsewander et alJJ2009l) . and perhaps most impor tantly, their 
locations relative to the host galaxies dFong et al.l 2010). As 
of mid-2010, X-ray and optical afterglows have been de- 
tected from 40 and 20 short GRBs, respectively, with the lat- 
ter sample providing accurate sub-arcsecond positions. Of 
these 20 events, 15 directly coincide with host galaxies with 
a wide distribution of appa rent magnitudes, and redshifts of 
zw 0.2-1 or b eyond (e.g., iBerger et alj[2007t lBergerll2009t 
iD'Avanzo et al.ll2009t) . However, the remaining five events 2 
do not appear to coincide with galaxies, and therefore provide 
an opportunity to assess the possibility of large progenitor off- 
sets, and to test the validity of the NS-NS merger models. 

Significant offsets have been claimed previously, in partic- 
ular for GRBs 050509b and 060502b wi th projected offsets 
of 39 ± 13 and 73 ± 19 kpc, respectively (iBloom et al.ll2006l 
2007). However, in both cases only X-ray positions are avail- 
able (3.5" and 4.4" radius, respectively), and the error circles 
contain sev eral galaxies consistent with a negligible offset 
(Berger et al. 2007; Bloom et al. 2007). Moreover, in the case 
of GRB 050509b the X-ray error circle intersected the outer 
regions of the host, raising the possibility that the progeni- 
tor system was formed in, rather than kicked to, the outskirts 
of the host. This possibility raises a crucial point, namely 
that a substantial physical offset from the center of the host 
does not necessarily point to a progenitor kick if th e burst still 
close ly coincides with the host light distribution (Fon g et all 
120101) . An illustrative example of this point is GRB 071227 
whose optical afterglow position coincides with the outskirts 
of edge-on disk galaxies, with an offsets of about 15 kpc from 
the host center (iD'Avanzo et a l. 2009; Fon g et al.ll2010l) . 

Large offsets have also been speculated in a few cases with 
precise optical afterglow positions (GRBs 06 1 201, 70809 , 
and 080503; IStratta et all 120071: iPerlev et al.l 120081 l2009h . 
However, these claims have not been investigated systemat- 
ically, mainly because they were treated on a case-by-case 
basis, with probabilistic arguments that prevented conclusive 
associations. These cases, combined with the ambiguity in- 
herent to X-ray positions, demonstrate that bursts with optical 
afterglows are essential for reaching any robust conclusions 
about progenitor offsets (due to kicks and/or a globular clus- 
ter origin). 

Here we present the first systematic study of short GRBs 
with optical afterglows and no coincident hosts, which com- 
bines their afterglow properties with the large-scale environ- 
ments. The purpose of this study is to statistically assess the 
possibility of offsets and to compare this with alternative ex- 
planations (e.g., a high redshift origin). To achieve this goal 
we set our study in the broader context of short GRBs that 
have optical afterglows and coincident hosts, as well as short 
GRBs with only X-ray positions. As we demonstrate through- 
out the paper, such a combined study is essential since offsets 
or high redshifts are expected to jointly affect both the after- 
glow properties and the large-scale environments. Our study 
also provides a robust statistical assessment of a posteriori 
chance coincidence probabilities, and the expected number of 
spurious associations, which cannot be properly assessed for 
individual bursts. 

2 These are GRBs 06 1201: IBergeil <2006j); IStratta et alJ J2007I); 
IFonget all 12010): 070809: Perlev et al. (2007, 2008); 080503~ |PerIev et alJ 
12009); 090305: Cenko et al. 12009); Berger & Kelson 12009); and 090515: 
Rowlinson etaLl 120101) 1 



The paper is organized as follows. In §|2] we present deep 
optical observations of the environments of GRBs 061201, 
070809, 080503, 090305, and 090515, as well as spectro- 
scopic observations of bright galaxies near the burst positions 
for GRBs 061201, 070809, and 090515. We study trends in 
the afterglow and prompt emission properties of bursts with 
and without coincident hosts in §|3] and determine a posteri- 
ori probabilities of chance coincidence as a function of pro- 
jected angular offset for galaxies near the position of each 
burst in S|4] We also determine projected physical and host- 
normalized offsets, and use these in conjunction with the af- 
terglow and prompt emission properties to address two sce- 
narios for the short bursts with optical emission and no coinci- 
dent bright hosts: (i) an origin in faint galaxies at z > 2 ( $4.21 ), 
or (ii) substantial offsets from galaxies at z ~ 0. 1 -0.5 or z > 1 
( j34.3l >. Finally, in §0 we present the offset distribution of all 
short GRBs with optical afterglows (in the context of scenario 
ii), and compare this distribution with predictions for NS-NS 
kicks and dynamically-formed NS-NS binaries. We further 
investigate whether the circumburst densities that are required 
for the measured optical magnitudes can be accommodated 
with a halo or IGM origin. We draw conclusions about the 
progenitors of short GRBs from our systematic study in S|6] 

2. SHORT GRB SAMPLE AND OBSERVATIONS 

We include in this investigation all 20 short GRBs with op- 
tical afterglow detections as of June 2010. This is the full 
subset of events for which sub-arcsecond positions are avail- 
able. We stress that this sample represents only about 1/3 of 
all short GRBs discovered to date, and about 1/2 of the sam- 
ple with X-ray afterglow detections. Thus, it is not a complete 
sample of short GRBs, but it is the only subset for which we 
can investigate the possibility of large offsets with meaningful 
statistical significance. As we demonstrate in §|5] we do not 
expect this sample to be strongly biased with respect to cir- 
cumburst density, at least for n > 10~ 5 cm" 3 . The properties 
of the 20 short GRBs, as well as events with deep optical lim- 
its, are summarized in Table Q] As can be inferred from the 
Table, some of the 20 events with only X-ray detections do 
not have optical follow-up observations, suggesting that the 
sample with optical afterglows may be largely representative. 

For the purpose of our investigation we define three sub- 
samples that will be used throughout the paper: (i) Sample 
1: short GRBs with detected afterglows and coincident host 
galaxies (15 bursts); (ii) Sample 2: the 5 short bursts with de- 
tected optical afterglows and no bright coincident hosts; and 
(iii) Sample 3: short GRBs with detected X-ray afterglows 
(from the Swift X-ray Telescope) but no optical detections de- 
spite rapid follow-up observations (11 bursts). 

2.1. Optical Imaging 

For the bursts in Sample 2 we use deep space- and ground- 
based optical observations to place limits on the brightness of 
underlying galaxies, and to assess the probability of chance 
coincidence for nearby galaxies. For GRB 061201 we use 
Hubble Space Telescope (HST) observations obtained with 
the Advanced Camera for Surveys (ACS) in the F814W fil - 
ter, with a total exposure time of 2244 s (Fong et~aT1l2010l) . 
For GRB 070809 we use r-band observations obtained on 
2008 January 14 UT with the Low Dispersion Survey Spec- 
trograph (LDSS3) mounted on the Magellan/Clay 6.5-m tele- 
scope, with a total exposure time of 1500 s. For GRB 080503 
we use HST Wide-Field Planetary Camera 2 (WFPC2) ob- 
servations in the F606W filter from 2009 January 30 UT, 



with a total exposure time of 4000 s dPerlev et all l2009h . 
We also use the limits on a coinci dent host inferred f rom a 
deeper stack of HST observations by Per lev et al.l(l2009l) . with 
m AB (F606W) > 28.5 mag. For GRB 090305 we use LDSS3 r- 
band observations obtained on 2010 May 8 UT with a total ex- 
posure time of 2400 s. Finally, forGRB 090515 we use r-band 
observations obtained with the Gemini Multi-Object Spectro- 
graph (GMOS) mounted on the Gemini-North 8-m telescope 
from 2009 May 15 UT with a total exposure time of 1800 s. 

The ground-based observations were reduced and analyzed 
using standard routines in IRAF. The ana lysis of the HST ob- 
servations is detailed in dFongetal .112010 1). The limiting mag- 
nitudes for all five observations are listed in Table |2] and im- 
ages of the five fields are shown in Figures [JJ-[5] 

2.2. Optical Spectroscopy 

In addition to the imaging observations, we obtained spec- 
troscopic observations of galaxies near the positions of GRBs 
061201, 070809, and 090515. Our observations of a galaxy 
located 16.3" from the afterglow position of GRB 061201 
(marked " S4" in Figu re [TJ| revealed a star forming galaxy at 
z = 0.111 dStratta et a l. 2007; Fong et al.|[2010h : see Figure[6] 

For GRB 070809 we obtained spectra of two galaxies lo- 
cated 5.9" and 6.0" from the optical afterglow position 
(marked "S2" and "S3", respectively in Figure O using 
LDSS3 on 2008 January 14 UT. The galaxy at a separation 
of 5.9" was previously identifie d as a star forming galaxy at 
z = 0.218 by lPerley etaD (l2008h . Here we find that the object 
at a separation of 6.0" is an early-type galaxy at z = 0.473, 
with no evidence for on-going star formation activity (Fig- 
ure |7). 

Finally, for GRB 0905 15 we obtained multi-object spec- 
troscopic observations with LDSS3 for nearly 100 galaxies 
within a 5' x 5' field centered on the GRB position. These ob- 
servations provide redshifts for several galaxies near the host 
position, including a star forming galaxy at z = 0.626 (5.8" 
offset; "SI" in Figure |5), an early-type galaxy at z = 0.403 
(14.0" offset; "S5"), and a star forming galaxy at z = 0.657 
(14.9" offset; "S6"); see Figure [8] In an upcoming paper we 
will demonstrate that the galaxy at z = 0.403 is a member of a 
cluster. 

3. AFTERGLOW PROPERTIES 

We begin our investigation by assessing the distribution 
of optical afterglow magnitudes for the three samples de- 
fined in ^2] The observed magnitudes and limits as a func- 
tion of time after the burst are shown in Figure [9] The me- 
dian observation time for the sample is about 8.5 hr after the 
burst. The distribution of detected optical afterglow magni- 
tudes is broad, ranging from tab ~ 21 to « 26 mag. The lim- 
its range from tab ^ 23 to > 25 mag (with the exception of 
GRB 080702 which on ly has a shallow limit of about 21 AB 
mag: Greco et al.ll2008h . 

The mean brightness and standard deviation for Sample 1 
are 3 (rp&) = 23.0 ± 1.3 mag. This is substantially fainter, 
by about an order of magnitude, than the afterglows of long 
GRBs on a comparable timescale (e.g., iKann et al] I2007L 
2008). It is also remarkably similar to the prediction of R rs 23 
mag at t ps 10 hr by (Pan aitescu et alJl2001[) . Since the avail- 
able limits are at least > 23 mag, we conclude that the bursts 

3 These numbers remain essentially unchanged if we extrapolate all mea- 
sured magnitudes to the fiducial time of 8.5 hr with a typical afterglow decay 
index of a = —0.75. 



lacking optical detections are drawn from a population with 
fainter afterglows. For the 5 bursts in Sample 2 we find a me- 
dian and standard deviation of (rAB) = 24.4 ± 1.4 mag, about 
1.4 mag fainter than the bursts with coincident hosts. 

The cumulative distributions for the three samples normal- 
ized to the fiducial time of 8.5 hr after the burst are shown in 
Figure [10] A Kolmogorov-Smirnov (K-S) test indicates that 
there is only an 8% probability that Sample 1 and Sample 2 
are drawn from the same underlying distribution of optical af- 
terglow brightnesses. Similarly, the probability that Sample 

1 and Sample 3 are drawn from the same underlying distri- 
bution is < 5% (an upper limit since the bursts in Sample 3 
are not detected in the optical). On the other hand, the prob- 
ability that Sample 2 and Sample 3 are drawn from the same 
underlying distribution is high, « 50%. 

The overall faintness of the optical afterglows in Sample 2 
and Sample 3 can be explained in two primary ways. First, 
they could result from systematically lower circumburst den- 
sities. In the standard afterglow model 4 with v m < v ovt < v c , 
the afterglow flu x scales as F v oc n l l 2 for a uniform medium 
(ISari et al.ll 998). Thus, adifferenceofabout+1.4magcanbe 
explained with a density that is about an order of magnitude 
lower than for the bursts in Sample 1. 

Alternatively, the fainter fluxes may be due to higher red- 
shifts for Sample 2 and Sample 3 compared to Sample 1 since 
the optical flux also scales as F v oc (1 + z) < - 3+p ^ 4 dj 2 , where di 
is the luminosity distance. The +1.4 mag difference corre- 
sponds to Az k, +0.5 (+1) for a Sample 1 mean redshift of 
z = 0.5 (z= 1). Similarly, the flux also depends on the total 
energy, with F v oc E < - 3+p ' > > 4 , and therefore a difference of +1 .4 
mag can be explained with an energy release lower by about 
a factor of 3. 

The various scenarios (lower density, lower energy, or 
higher redshift) can be further explored through a compari- 
son of the prompt 7-ray emission, and the relation between 
the optical and X-ray afterglow brightness. In the framework 
of the standard GRB model we do not expect lower densities 
to impact the prompt emission since it is expected to be pro- 
duced by internal processes (shocks or magnetic dissipation) 
that do not depend on the external medium. On the other hand, 
lower energies or higher redshifts will tend to systematically 
affect the prompt and afterglow emission. 

In Figure [TU we plot the distributions of 7-ray fluence (F 7 ), 
afterglow X-ray flux at the fiducial time of 8 hr (Fx$), and 
duration (T90) for the three samples. These distributions allow 
us to explore the underlying reason for the difference in opti- 
cal afterglow brightness between the three samples. First, we 
find that the distribution of F 1 values for Sample 1 has a mean 
value that is about a factor of 5 times larger than for Sample 

2 and Sample 3. This is indicative of higher redshifts for the 
latter two samples if the isotropic-equivalent energies of all 
short GRBs are similar, or alternatively a lower energy scale 
if the redshifts are similar. 

In the same vein, we find that the distributions of Fx.% val- 
ues for Sample 2 and Sample 3 have lower means than for 
Sample 1. This result is indicative of overall fainter afterglow 
emission for the former two samples, and this can again be 
explained in the context of lower energies or higher redshifts. 

4 We use the standard synchrotron spectrum definitions: v m is the charac- 
teristic synchrotron frequency corresponding to electrons with the minimum 
Lorentz factor (7,,,) of the electron distribution, N(j) oc j~P; p us ually has a 
value of ss 2.2 — 2.5; and u c is the synchrotron cooling frequency jSari et all 



Unlike in the case of the 7-ray emission, a lower circum- 
burst density would also account for the fainter X-ray fluxes 
if vx < v c . Finally, we find that the durations of the bursts 
in Sample 2 and Sample 3 are shorter by about a factor of 2 
compared to the events with optical afterglows and coincident 
hosts, although there is substantial scatter in all three samples. 
The shorter durations are not trivially explained in the context 
of lower energies, higher redshifts, or lower densities. 

To explore this result in more detail we plot the observed 
fiuence as a function of duration for the events in all three 
samples (Figure [121 ). We find that there is a mild positive 
correlation between the two quantities, but that the events in 
Sample 1 appear to have larger fiuences at a given duration 
compared to the events in Sample 2 and Sample 3. This is 
indicative of lower 7-ray fluxes for the latter two samples, 
possibly as a result of higher redshifts. Higher redshifts will 
also shift the intrinsic durations of the bursts in Sample 2 and 
Sample 3 into better agreement with the bursts in Sample 1. 

We therefore conclude that the differences in prompt emis- 
sion and optical/X-ray afterglow properties are consistent 
with a higher redshift origin for the bursts in Sample 2 and 
Sample 3. The fainter afterglow emission is also consistent 
with lower density environments for these two samples, al- 
though this does not clearly explain the differences in prompt 
emission (at least in the framework of the standard GRB 
model). We return to the discussion of low density versus a 
high redshift origin in $4]when we investigate the host galaxy 
properties. 

In addition to the overall faintness of the optical and X-ray 
afterglows, a substantial difference in density may also be im- 
printed on the ratio of optical to X-ray brightness. This is be- 
cause the synchrotron cooling frequency depends on density 
as v c oc ri~ l , and is therefore expected to transition across the 
X-ray band as the density decreases. For v c > vx the X-ray 
and optical bands occupy the same portion of the synchrotron 
spectrum, with a resulting spectral index of j3 = -{p- l)/2 « 
-0.6 to fh -0.75 while for v c < vx (i.e., high density), the 
spectrum between the two bands will be steeper, reaching a 
maximum value of m -1.25 when v c w v opt . In Figure Q~3] 
we plot the X-ray flux versus optical magnitude for all three 
samples. For each burst the fluxes in the optical and X-rays 
are taken at similar times after the burst, or extrapolated to 
a common time. The correction factors due these extrapola- 
tions are generally 5 < 2 (Table [TJ. For the combined Sample 
1 and Sample 2 we find a clear correlation between the fluxes 
in the two bands, leading to a mean optical to X-ray spectral 
index of (/3ox) = -0.72 ±0.17. This is essentially indistin- 
guishable from the ra tio for long GRBs, (/3ox) = -0.65 ±0.35 
( Jakob sson et alj 120041) . The median values are consistent 
with v c > vx, but exhibit dispersion that is likely due to scat- 
ter in the values of p and/or the location of v c relative to the 
X-ray band. 

The similarity of /?ox for long and short GRBs does not 
necessarily indicate that the densities are similar for the two 
samples. In particular, if v c is located close to the X-ray 
band for long GRBs, while for short GRBs v c ^> vx (due to a 
lower density), the effect on /?ox will be marginal, particularly 
within the overall observed scatter. For example, with p = 2.5 
the difference in /3ox between a model with v c exactly inter- 
mediate between the optical and X-ray bands, and a model 
with v c > vx, is A/?ox ~ 0.25. On the other hand, the scatter 
resulting from a range of p = 2.2-2.5 is of the same order, 

5 We do not extrapolate X-ray upper limits. 



A/3ox = 0.15. Similarly, the nearly equivalent median /3ox 
values may indicate that for both GRB populations v c > vx- 
In this case, the resulting lower limits on v c therefore prevent 
the use of /?ox as an indicator of density. 

Comparing Sample 1 and Sample 2, we find no clear dif- 
ference in Pox (Figure [T3l. The same is true for the bursts 
in Sample 3, which are all consistent with the same relation 
given the optical upper limits and a mix of X-ray detections 
and upper limits. Thus, the ratio of optical to X-ray flux does 
not allow us to distinguish redshift/energy and density effects 
between the three samples. 

To summarize, the optical afterglows of short GRB without 
coincident hosts (or with only optical limits) are systemati- 
cally fainter than those of short GRBs with coincident hosts. 
The same is true for their X-ray fluxes and 7-ray fiuences. 
The fainter afterglows may reflect lower densities (by an order 
of magnitude), but this does not naturally explain the lower 
7-ray fiuences. Alternatively, the fainter afterglows and 7- 
ray fiuences can be explained as a result of higher redshifts 
(Az ~ 0.5 - 1) or lower energies (by about a factor of 3). 

4. LARGE-SCALE ENVIRONMENTS 

We next turn to an analysis of the large-scale environments 
of the bursts in Sample 2, partly in comparison to the hosts 
of bursts in Sample 1. As indicated in ^2) we place limits of 
tab ~ 25.4-26.5 mag on the brightness of any galaxy under- 
lying the five short GRB positions (Tabled. These limits are 
several magnitudes fainter than the measured brightnesses of 
short GRB hosts at z < 1 (Sample 1). 

4.1. Probabilities of Chance Coincidence 

To assess the potential that galaxies near each of the 5 bursts 
are the hosts, we calculate their probability o f chance coinci- 
dence. We follow the methodology of Blo om et al.l {2002), 
namely, we determine the expected number density of galax- 
ies brighter than a measured magnitude, m, using the results of 
deep optical galaxy surveys (Hogg et al. 1997l lBeckwifh et all 
2006): 
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0.33xln(10) 

The probability for a given separation, P(< 8R), is then given 
by 

P«^)=l-e- 7r(5R)2<T( ^" ) , (2) 

where we use the fact that for offsets substantially larger than 
the galaxy size, 8R is the appropriate radius in Equation [2] 
(iBloom et al.ll2002l) . 

The resulting distributions for each field are shown in Fig- 
ure [14] We include all galaxies that have probabilities of 
< 0.95. We find that for 4 of the 5 bursts, faint galaxies 
(~ 25-26 mag) can be identified within « 1.6-2" of the 
afterglow positions, with associated chance coincidence prob- 
abilities of « 0.1 -0.2; in the case of GRB 090515 we do not 
detect any such faint galaxies within « 5" of the afterglow 
position. For GRB 080503 we also include the galaxy at an 
off set of 0.8" and m A B (F606W) = 27.3 ±0.2 mag identified 
by iPerley et aD, (2009) based on their deeper stack of HST 
observations. On the other hand, for 4 of the 5 bursts we 
find that the galaxies with the lowest probability of chance 
coincidence, rs 0.03-0.15, are brighter objects with offsets 
of about 6-16" from the burst positions; only in the case of 
GRB 080503 the lowes t chance coinciden ce is associated with 
the nearest galaxy (see lPerlev et al.l 2009). 



For comparison we repeat the same analysis for short 
GRBs from Sample 1 which have faint coincident hosts 
(GRB 060121: 26.0 AB mag; GRB 060313: 26.6 AB mag; 
and GRB 070707: 27.3 AB mag). The results of the prob- 
ability analysis are shown in Figure Q3] We find that in all 
three cases, the coincident hosts exhibit the lowest probabil- 
ity of chance coincidence, « 0.02-0.05. Only in the case 
of GRB 070707 do we find galaxies with SR > few arcsec 
that have P(< 5R) < 0.1. Thus, these three bursts are con- 
sistent with negligible offsets from faint galaxies, presumably 
atz> 1. 

The use of a posteriori probabilities to assign unique galaxy 
associations is fraught with difficulties. First, for a given 
apparent brightness, galaxies located further away from the 
GRB position, potentially due to larger kicks and/or longer 
merger timescales in the NS-NS merger framework, have 
higher probabilities of chance coincidence. Since we have no 
a priori model-independent knowledge of the range of possi- 
ble kicks and merger timescales, we cannot rule out galaxies 
at very large offsets for which P{< 5R) ~ 1. Indeed, a rea- 
sonable constraint of vyck ^ 10 3 km s" 1 and T mergel - < 10 Gyr 
leads to only a weak constraint on the offset of < 10 Mpc. At 
z = 0.1 (z = 1) this corresponds to about 1.5° (0.3°), a pro- 
jected distance at which nearly all galaxies will have a chance 
coincidence probability of order unity. 

A second difficulty is that we are using angular offsets, 
which ignore the potential wide range of redshifts (and by 
extension also luminosities) of the various galaxies. For ex- 
ample, if the faint galaxies with small offsets are located at 
z > 1, the corresponding physical offsets are ~ 15 kpc, while 
if the galaxies at ~ 10" offsets are located at z ~ 0.3 (see ©, 
the offsets are only somewhat larger, ~ 30 kpc. A galaxy at 
an even lower redshift, z ~ 0.1, with an offset of 50 kpc will 
be located about 30" from the GRB position and incur a large 
penalty in terms of chance coincidence probability. It is im- 
portant to note, however, that galaxies at lower redshift will 
generally have brighter apparent magnitudes, partially com- 
pensating for the larger angular separations (Equations Q] and 
|2). In only a single case (GRB 070809) we find a galaxy with 
P(< SR) < 0.1 at 6R > 1' (which at z = 0.043 for this galaxy 
corresponds to a physical offset of about 100 kpc). 

A final complication, which is not unique to this subset of 
events, is that we can only measure projected offsets, 5R = 
SRj,u x cos(0). The measured offsets can be used as lower 
limits on the actual offsets, while for the overall distribution 
we can apply an average correction factor of tt/2, based on 
the expectation value for the projection factor, cos(#). 

Despite these caveats we can address the probability that all 
of the associations are spurious. This joint proba bility is sim- 
ply th e product of the individual probabilities ( Bloom et al.1 
l2002|) . For the faint galaxies at small angular separations the 
probability that all are spurious associations is P a n ss 8 x 10~ 5 , 
while for the galaxies with the lowest probability of chance 
coincidence the joint probability is nearly 30 times lower, 
P a u w3x 10~ 6 . Conversely, the probabilities that none of 
the associations are spurious are k, 0.42 and w 0.59, respec- 
tively. These values indicate the some spurious coincidences 
may be present for Sample 2. Indeed, the probabilities that 
1, 2, or 3 associations are spurious are [0.40, 0.15, 0.027] and 
[0.34, 0.068, 0.006], respectively. These results indicate that 
for the faint galaxies it is not unlikely that 2-3 associations 
(out of 5) are spurious, while for the brighter galaxies 1-2 
associations may be spurious. This analysis clearly demon- 
strates why a joint statistical study is superior to case-by-case 



attempts to associate short GRBs with galaxies at substantial 
offsets. 

We therefore conclude that despite the weaknesses inherent 
to a posteriori probabilities, we conclude that there is stronger 
statistical support for an association of at least some of the 
5 bursts in Sample 2 with bright galaxies at separations of 
~ 10", than for an association with the faint galaxies at sepa- 
rations of ~ 2" . Clearly, we cannot rule out the possibility that 
in reality the hosts are a mix of faint and bright galaxies with a 
range of angular offsets of > 2". We note that if deeper obser- 
vations eventually lead to the detection of underlying galax- 
ies (< 0.5") at the level of « 27 mag, the associated chance 
coincidence probabilities will be k, 0.05 per object, and the 
joint probabilities will be only slightly higher than for the 
bright galaxies with ~ 10" offsets. On the other hand, if we 
can achieve magnitude limits of > 28 mag on any coincident 
hosts, the resulting probabilities of chance coincidence will 
be larger than for the offset bright galaxies. Thus, eliminating 
the possibility of underlying hosts at the level of > 28 mag 
is of the utmost importance. So far, only GRB 080503 has 
been observed to such a depth with no detected coincident 
host (iPerlev et al.ll2.Q09b . but observations of the full sample 
are requ ired for a robust statistical comparison. As we discuss 
in {34.2l below. such deep limits will also reduce the probabil- 
ity of underlying hosts based on redshift arguments. 

Beyond the use of projected angular offsets, we note that 
the faint galaxies near the burst positions are likely to have 
projected physical offsets of about 6 15 kpc. To assess the pro- 
jected physical offsets for the galaxies with the lowest proba- 
bility of chance coincidence we measured spectroscopic red- 
shifts in three cases (GRBs 061201, 070809, and 090515; ©. 
In the case of GRB 061201, the galaxy redshift o f z = 0-111 
leads to a projected physical offset of 32.4 kpc (Fo nget al.l 
120101) . In the case of GRB 070809, the lowest probability 
of chance coincidence is not associated with the star form- 
ing galaxy 7 at z = 0.218 identified by Perley et al. (2008), but 
instead belongs to the early-type galaxy identified here, which 
has a redshift of z = 0.473, and hence an offset of 34.8 kpc. 
Finally, for GRB 090515, the lowest probability of chance 
coincidence is associated with a cluster early-type galaxy at 
z = 0.403, leading to a physical offset of 75 kpc. 

Thus, two scenarios emerge from our investigation of the 
large-scale environments: (i) the bursts are spatially coinci- 
dent with currently-undetected galaxies (with tab ^ 26 mag); 
or (ii) the bursts have substantial offsets of at least ps 15-75 
kpc depending on whether they are associated with faint 
galaxies at small angular separations, or brighter galaxies at 
Z ~ 0.1 -0.5; for the ensemble of 5 events the larger offsets 
are statistically more likely than the ~ 15 kpc offsets. In the 
context of large offsets, even larger values may be possible if 
the bursts originated in galaxies with larger separations and 
P(<5R)~\. 

These two scenarios echo the possibilities that emerged 
from the analysis of the afterglow and prompt emission prop- 
erties (©. Distinguishing between these possibilities is 
clearly of fundamental importance to our understanding of 
short GRBs: the former scenario will point to a population of 
very faint hosts (likely at high redshifts), while the latter sce- 
nario will provide evidence for large offsets (due to kicks or a 

6 Although the redshifts of these galaxies are not known, the angular di- 
ameter distance is nearly independent of redshift beyond z ~ 1 , which is ap- 
propriate for these faint host galaxies. 

7 We note that even if the burst was associated with this galaxy, the corre- 
sponding offset would be 20.6 kpc. 



globular cluster origin) and hence NS-NS/NS-BH progenitors 
for at least some short GRBs. 

4.2. Scenario 1: Undetected Faint Hosts at High Redshift 

We can place upper limits on the redshifts of the GRBs 
in Sample 2 based on their detections in the optical (i.e., 
the lack of complete suppression by the Lya forest). The 
afterglow of GRB 061201 was detected in the ultraviolet 
by the Swift/UVOT and it is therefore located at z < 1.7 
dRoming et alj 12006). The remaining four bursts were de- 
tected in the optical g- or r-band, and can therefore be placed 
at z < 4 (Tabled}- 

We place additional constraints on the redshifts of any 
underlying hosts using the existing sample short GRB host 
galaxies. In Figure [16] we plot the r-band magnitudes as a 
function of redshift for the all available short GRB hosts from 
Sample 1. For the faint hosts without known redshifts (GRBs 
060121, 060313, and 070707), we place upper limits on the 
redshift using optical detections of the afterglows (Table [TJ. 
A wide range of host magnitudes, tab ~ 16.5-27.5 mag, is 
appar ent. We also plot the r-band magnitudes of long GRB 
hosts ( JSavaglio et alj|2009l) . as well as the r— z phase space 
that is traced by galaxies with luminosities of L = 0. 1 - 1 L* . 
We use the appropriate value of L* as a function of redshift, 
taking into account the evolving galaxy luminosity function 
(Steidel et al. 1999; Blanton et al. 2003; Wi llmer et al] 2006; 
Reddv & Steidel 2009). We find excellent correspondence be- 
tween the hosts of long and short GRBs, and the phase-space 
traced by 0. 1 - 1 L* galaxies, at least to z ~ 4. In the context 
of these distributions, the available limits for the short GRBs 
in Sample 2 translate to redshifts of z > 1 .5 if they are 0. 1 L* 
galaxies, or z > 3 if they are L* galaxies. The latter lower lim- 
its are comparable to the redshift upper limits inferred from 
the afterglow detectio ns. We note that fo r GRB 080503, the 
limit of > 28.5 mag dPerlev et alj 2009) places even more 
stringent limits on the redshift of an underlying 0.1-1 L* 
galaxy. 

The redshifts of z > 1.5 for putative 0.1 L* hosts are 
consistent with the faintness of the optical afterglows, from 
which we inferred Az sw 0.5 - 1 compared to Sample 1 (©. 
We note, howe ver, that the one kn own short GRB at z > 2 
(GRB 090426; lLevesqueet~aT]|2010l) has a host galaxy lumi- 
nosity of ~ 2 L*, which may suggest that the appropriate red- 
shift lower limits are z > 3. 

The possibility that the five bursts originated at z > 3 leads 
to a bimodal redshift distribution (Figure [TBI). Nearly all of 
the bursts in Sample 1 with a known redshift (9/10) have 
Z p» 0.2- 1, with a median of (z) ~ 0.5; the sole exception 
is GRB 090426 at z = 2.61. The three bursts with faint coin- 
cident hosts have upper limits of z < 4 from afterglow detec- 
tions, while lower limit of z > 1.5-2 can be placed on these 
hosts if they have L > 0. 1 L* . Adding the Sample 2 bursts 
with the assumption that they have z > 3 will furthermore re- 
sult in a population of short GRBs with a median of z ~ 3, and 
leave a substantial gap at z ~ 1-2 (FigurefToTl. If the 5 bursts 
are instead hosted by 0. 1 L* galaxies, the inferred lower limits 
on the redshifts (z > 1 .5) lead to a potentially more uniform 
redshift distribution. 

It is difficult to explain a bimodal redshift distribution with a 
single progenitor population such as NS-NS binaries, without 
appealing to, for example, a bimodal distribution of merger 
timescales. Another possibility is two distinct progenitor pop- 
ulations, producing bursts of similar observed properties but 
with distinct redshift ranges. While these possibilities are dif- 



ficult to exclude, they do not provide a natural explanation for 
the short GRB population. 

A final alternative explanation is that any underlying hosts 
reside at similar redshifts to the known hosts in Sample 1 (z ~ 
0.5), but have significantly lower luminosities of < 0.01 L*. 
This scenario would not naturally explain why the bursts in 
Sample 2 have fainter optical and X-ray afterglows, as well 
as lower 7-ray fluences. We therefore do not consider this 
possibility to be the likely explanation. 

4.3. Scenario 2: Large Offsets 

While higher redshifts may explain the lack of detected 
hosts, the fainter afterglows, and the weaker 7-ray fluences 
of the bursts in Sample 2, this scenario suffers from several 
difficulties outlined above. The alternative explanation is that 
the bursts occurred at significant offsets relative to their hosts, 
and hence in lower density environments that would explain 
the faint afterglow emission (though possibly not the lower 7- 
ray fluences). As we demonstrated in 34.11 the offsets may be 
~ 2" (~ 15 kpc) if the bursts originated in the faint galaxies at 
the smallest angular separations, or ~ 10" (~ 30-75 kpc) if 
they originated in the brighter galaxies with the lowest prob- 
ability of chance coincidence. Below we address the implica- 
tions of these two possible offset groups through a compari- 
son to the offsets measured for the bursts in Sample 1 (e.g., 
iFong et al.ll20Toh . 

We plot the distributions of projected angular offsets for the 
short GRBs with and without coincident hosts in Figure [T7l 
The offsets for Sample 1 have a mean and standard deviation 
of about 0.7 ± 0.7" and a range of about 0. 1-3". Modeled 
with a log-normal distribution, the resulting mean and width 
in units of arcseconds are log((5/?) f=a -0.2 and <j\ 0% {$r) ~ 0.35. 
If we associate the short bursts in Sample 2 with the galax- 
ies that have the lowest probability of chance coincidence, 
the resulting distribution has (SR) = 9 ± 6 arcsec. This is 
clearly distinct from the distribution for Sample 1. The ef- 
fect is less pronounced if we associate the bursts with the 
galaxies located at the smallest angular offsets (Figure \YT\ . 
Even for these separations the mean and standard deviation 
are \og(8R) w 2.6 ± 1 .8 arcsec. 

As noted in 34.11 the projected angular distances may not 
be the most robust quantity for measuring the offset distribu- 
tion. An alternative qu antity is the offse t normalized by each 
host's effective radius dFonget alj|20 10). This quantity takes 
into account the varying sizes of the hosts due to both intrin- 
sic size variations and redshift effects. It also gives a better 
indication of whether the burst coincides with the host light 
or is significantly offset. As shown in Figure [18] the host- 
normalized offsets of Sample 1 have a mean and standard de- 
viation of about 1 ±0.6 R e , and a range of about 0.2-2^. A 
log-normal fit results in a mean of log{SR/R e ) « and a width 
°f ciogC&R/R,) ~ 0-2- The bursts in Sample 2 have much larger 
host-normalized offsets, with (8R/R e ) = 7.3 ±2.3 if they origi- 
nated in the galaxies with the lowest chance coincidence prob- 
ability. Even if we associate the bursts with the nearest faint 
hosts, the distribution has a mean of about 4 R e , reflecting the 
fact that the effective radii of the faint galaxies are smaller 
than those of the brighter ones. 

Finally, we plot the projected physical offsets in Figure [T9l 
The mean and standard deviation for Sample 1 are SR = 4.2 ± 
3.8 kpc, and a log-normal fit results in a mean of log(^) « 0.5 
and a width of <r\ og (SR) « 0.3. On the other hand, the bursts in 
Sample 2 have a mean offset of about 19 kpc if they arise in 
the faint galaxies with small angular separation, or about 40 



kpc if they arise in the brighter galaxies, again pointing to 
distinct distributions. 

The distributions of angular, physical, and host-normalized 
offsets exhibit a clear bimodality if we associate the bursts in 
Sample 2 with the galaxies at z ~ 0. 1-0.5. This is particularly 
apparent in the more meaningful quantities, namely physical 
and host-normalized offsets (Figures [T8l and |T9b. The effect 
is still apparent, though less pronounced in the case of asso- 
ciation with the faint galaxies at z > 1. Thus, if the offset 
scenario is correct, the resulting distributions point to a possi- 
ble bimodality rather than a single continuous distribution of 
offsets. 

The cumulative distributions of physical offsets for Sample 
1 alone, and in conjunction with the two possible offset groups 
for Sample 2 are shown in Figure[20] The combined distribu- 
tions have a median of about 4 kpc, driven by the bursts with 
coincident hosts. However, there is a clear extension to larger 
physical offsets in the case of association with the brighter 
galaxies, with about 20% of all objects having SR > 30 kpc. 
The cumulative distributions are particularly useful for com- 
parison with NS-NS merger models since predictions exist for 
both the kick scenario and the globular cluster origin model. 
We turn to this discussion below. 

5. DISCUSSION AND IMPLICATIONS 

We have shown that the lack of host galaxy detections in 
coincidence with the bursts in Sample 2 indicates either a 
high redshift (z > 2) origin, or substantial projected offsets of 
~ 15-75 kpc (xi4-lOR e ). Both scenarios can account for the 
lack of underlying galaxy detections and the fainter optical/X- 
ray afterglows (either due to redshift or density effect); the 
possibility of large offsets from galaxies atz ~ 0.1 -0.5, how- 
ever, does not naturally explain the lower 7-ray fiuences and 
somewhat shorter durations of the bursts in Sample 2. Indeed, 
the lower fiuences can be more easily accommodated in the 
case of undetected high redshift coincident hosts, or in the 
case of faint hosts with offsets of ~ 15 kpc. 

Below we test the results of the offset scenario against pre- 
dictions for NS-NS merger models and dynamically-formed 
NS-NS binaries in globular clusters. We also investigate 
whether the densities expected at offsets of ~ 15-75 kpc (i.e., 
galaxy halos or the intergalactic medium) can accommodate 
the observed optical magnitudes. Finally, we discuss the im- 
plications of the various scenarios for short GRB energetics. 

5.1. Comparison to NS-NS Merger Models 

The sample of short GRBs with optical afterglows repre- 
sents about 1/3 of all short bursts, and may thus not be fully 
representative. One often-discussed bias is that the bursts 
with optical afterglows require a high circumburst density, 
and therefore have negligible offsets. However, from our anal- 
ysis in this paper it is clear that one explanation for the lack 
of coincident hosts for the bursts in Sample 2 is indeed large 
offsets, despite their detection in the optical band. 

The bursts with only X-ray afterglow detections 
(Swift/XRT) have typical positional uncertainties of 
ax ~ 2-5" and therefore lead to a deeper ambiguity 
about the identity of the hosts. We have previously argued 
that most of these bursts are associated with galaxies con- 
sistent with a negligible offse t (or a s large as tens of kpc; 
iBerger et~al . 2007; Fong et al] 120101) . In some cases larger 
offsets have been advocated based on a poste riori chance 
coinc idence probabil ities (e.g., GRB 05 0509b: Bloo m et alj 
12005 GRB 060502b: iBloom et alJl2007h . However, from the 



analysis in 34. H it is clear that while the inferred probabilities 
are only weakly dependent on the positional accuracy when 
SR > o~x (i.e., galaxies at large offsets), a large penalty is 
incurred for faint galaxies located within the X-ray error 
circles since in that case the appropriate value in Equation [2] 
is SR = ax- As a result, the faint galaxies will have chance 
coincidence probabilities of P(< SR) ~ 1, To avoid this 
complication in our comparison to NS-NS model predictions 
we restrict the analysis to the sample with optical afterglow 
positions. 

In recent work by Fon g et all J2010h we have shown that 
the offset distribution for short GRBs with optical and/or 
X-ray positions, and accounting for incompleteness due to 
bursts with only 7-ray positions, is broadl y consistent with 
the pre dictions of N S-NS merger models dBloom et alJI 19991 ; 
Frveretal. 1999; Belczynski et al. 2006). In particular, we 
concluded that > 25% of all short bursts have SR < 10 kpc, 
while > 5% have offsets of > 20 kpc. We repeat the anal- 
ysis here using the offsets inferred in 34.31 for both the faint 
galaxies at small angular separations and the brighter galaxies 
at larger separations. As shown in Figure |20l the model pre- 
dictions have a median of about 6 kpc, compared to about 4 
kpc for the observed sample. On the other hand, the models 
predict 10-20% of offsets to be > 30 kpc, in good agree- 
ment with the observed distribution in both the ~ 15 kpc and 
~ 40 kpc scenarios. We note that the overall smaller offsets 
measured from the data may be due to projections effects. In- 
deed, the mean correction factor of tt/2 nicely reconciles the 
theoretical and observed distributions. 

in the previous section we noted a bimodality in the phys- 
ical and host-normalized offsets for Sample 1 and Sample 
2 (Figures [18] and [19]>. In the framework of NS-NS bi- 
nary kicks this bimodality may indicate that the binaries gen- 
erally remain bound to their host galaxies, thereby spend- 
ing most of their time at the maximal distance defined by 
dmax = 2GMh os t/ vy^v (i.e., wit h their kinetic energy stored as 
potential energy; Bl oom et alj [2007). This would require typ- 
ical kick velocities of less than a few hundred km s" 1 . 

We further compare the observed offset distribution to pre- 
dictions for dynamically-formed NS-NS binaries in globu- 
lar clusters, with a range of host galaxy virial masses of 
5 x 10 10 - 10 12 M dSalvaterra et a l. 2010). These models 
predict a range of only w 5-40% of all NS-NS mergers to 
occur within 10 kpc of the host center, in contrast to the ob- 
served distribution with about 70% with SR < 10 kpc. We 
stress that this result is independent of what offsets we assign 
to the bursts in Sample 2 since they account for only 1/4 of the 
bursts with optical afterglows. On the other hand, the globu- 
lar cluster origin may account for the bimodality in the physi- 
cal and host-normalized offsets (Figures [T8land[T9l), with the 
objects in Sample 2 arising in globular clusters and the ob- 
jects with coincident hosts arising from primordial NS-NS bi- 
naries. This possibility also agrees with the predicted frac- 
tion of dynamically- formed NS-NS binaries of ~ 10-30% 
dGrindlav et al.l l2006). The cumulative offset distributions for 
Sample 2 alone (assuming the hosts are the galaxies with the 
lowest probability of chance coincidence) is well-matched by 
the range of predictions for dynamically-formed NS-NS bina- 
ries in globular clusters (Figurel20l>. In this scenario, however, 
the implication is that short GRBs outside of globular clusters 
do not experience kicks as expected for NS-NS binaries since 
the largest measured offset is only 15 kpc. 

Unless the populations of short GRBs with only X-ray or 
7-ray positions have fundamentally different offset distribu- 
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tions, we conclude that the measured offsets of short GRBs 
and the predicted offsets for NS-NS kicks are in good agree- 
ment, if we treat all short GRBs with optical afterglows as 
a single population. Alternatively, it is possible that the bi- 
modal distributions of physical and host-normalized offsets 
point to a progenitor bimodality, with the bursts in Sample 2 
originating in globular clusters. 

5.2. Circumburst Densities 

Large offsets are expected to result in low circumburst den- 
sities, and we now investigate whether the observed optical 
magnitudes for Sample 2 can be used to place meaningful 
constraints on the offsets. The optical afterglow magnitudes 
do not provide a direct measure of the circumburst density 
since they depend on a complex combination of density, en- 
ergy, and fractions of the burst energy imparted to the radiat- 
ing electrons (e e ) and magnetic fields (eg). From our compar- 
ison of the X-ray and optical afterglows (@ we were unable 
to clearly locate the cooling frequency, which in principle can 
provide additional constraints on the density. However, we 
can still gain some insight into the required circumburst den- 
sities using a few basic assumptions. We assume that the j- 
ray energy (£ 7 ; s0 ) is a reasonable proxy for the total energy, 
that e e ,€B < 1/3, and that p = 2.2. We further use a fiducial 
redshift of z = 0.5 (Figure [ToTl and a fiducial observation time 

of8hr(®- 

To determine £ 7] i S o we use the 15- 150 keV fiuences mea- 
sured by Swift (Table [TJ, and determine a correction factor 
to account for the incomplete energy coverage using bursts 
that have also been detected by satellites with a broader en- 
ergy range (Figure [121). We find a typical correction factor of 
w 5 for the « 10— 10 3 keV range. For the fiducial redshift of 
Z = 0.5, the mean energy for short GRBs with optical after- 
glows is is-y.iso w 2 x 10 51 erg, when we include the correction 
factor (Figure |2"TV 

With these assumptions we find the following relation be- 
tween the optical b rightness and the circumburst density 
dGranot & Sarill2002h : 



Fv, opt < 1 mJy x n 



1/2 



(3) 



Given the typical observed fluxes of about 2 /iJy for Sam- 
ple 1, and about 0.6 fiJy for Sample 2, we infer typical min- 
imum densities of > 4 x 10~ 6 and >4x 10~ 7 cm" 3 , respec- 
tively. Thus, the observed optical brightnesses can be pro- 
duced even at very low densities that are typical of the IGM. 
This indicates that we cannot rule out the large offset sce- 
nario based on density arguments. Indeed, even if we allow 
for more typical values of e e « e B « 0. 1 , the resulting densities 
are « 4 x 10~ 4 and 4 x 10~ 5 cm -3 , respectively. These results 
compare favorably with predictions for NS-NS kicks, which 
suggest that most mergers will occur at densitie s of > 10~ 6 
cm -3 (Perna & Belczvnski 2002; Belczvnski et al. 2006). The 
distribution of densities in globular clusters is not well known, 
preventing a meaningful compar ison to our inferred minimum 
densities (ISalvaterra et al.ll2010h . 

5.3. Energetics 

Finally, we investigate the energetics of the bursts in Sam- 
ple 2 in the context of a high-redshift origin (z > 3), an ori- 
gin in the faint galaxies at w 2" separations (z ~ 1), and an 
origin in the galaxies with the lowest probabilities of chance 
coincidence (z w 0.1 -0.5). The resulting values of E 7 j so 
are shown in Figure |2T] For the lowest redshift origin, the 



inferred values are ps (1-5) x 10 49 erg, for az~l origin 
they are «5x 10 49 - 10 51 erg, and for a z > 3 origin they are 
«5x 10 50 -5x 10 5I erg. 

For comparison, the mean value for the bursts in Sample 1 
is (£ 7l i S0 ) « 8 x 10 50 erg, with a range of « 10 49 -5 x 10 51 
erg. However, there is a clear redshift dependence for the 
measured E 1 ,\ so values (Figurel2T1i. with w5x 10 49 erg at z < 
0.5 and w 5 x 10 50 erg at z > 0.5. Thus, the bursts in Sample 
2 generally fit within the known range of isotropic energies 
regardless of their actual redshift. Indeed, at z ~ 3 they exhibit 
similar values of E 7 i s0 to that of GRB 090426. As a result, we 
cannot use the inferred energy release as a clear discriminant 
of the redshift range for Sample 2. 

If the bursts in Sample 2 indeed originated at high redshifts, 
the resulting isotropic-equivalent energies suggest that there 
is either a large spread in the energy release of short GRBs (at 
least 2 orders of magnitude) or a large variations in the ejecta 
geometry. If the typical energy for short GRBs is about 5 x 
10 49 erg, as indicated by the nearest events 8 , then the inferred 
values of E 7 j so w2x 10 51 at z ^ 3 indicate opening angles 
of m 10- 15°. This is similar to the op ening angle of about 
7° that was inferred for GRB 051221 (iBurrowset all 120061; 
ISoderberg et~alll2006l) . 

6. CONCLUSIONS 

We undertook the first systematic study of short GRBs with 
detected optical afterglows (and hence sub-arcsecond posi- 
tion) but no coincident host galaxies to limits of tab > 26 
mag. We find that the optical afterglows of these bursts are 
fainter by about 1.4 mag compared to the optical afterglows 
of short GRBs with coincident hosts. They are similarly 
fainter in X-rays, and have somewhat lower 7-ray fiuences 
and slightly shorter durations. Both samples have similar ra- 
tios of X-ray to optical flux, which are moreover similar to 
the ratios measured for long GRBs. The fainter afterglows of 
the bursts lacking coincident hosts may be due to lower den- 
sities, lower energies, or higher redshifts. However, we note 
that only the scenarios with lower energies or higher redshifts 
naturally explain the faintness of the prompt emission in the 
context of the standard GRB model. This is because in the 
context of internal processes (shocks or magnetic dissipation) 
the external density should not play a role. 

We further use deep optical imaging to determine the prob- 
ability of chance coincidence for galaxies in the field around 
each burst, and to place redshift limits for underlying hosts 
under the assumption that they are drawn from the same dis- 
tribution of the detected short GRB hosts (Lw0.1-1 L*). 
This analysis leads to the following possible scenarios: (i) the 
underlying hosts are fainter than ~ 26 mag, indicative of red- 
shifts of z > 1.5 (if L ~ 0.1 L*) or z > 3 (if L ~ L*); or (ii) 
the hosts are galaxies with substantial offsets — either faint 
galaxies at separations of w 2" (« 15 kpc for z > 1) or brighter 
galaxies, which we find to be located at z ~ 0.1 -0.5, with 
offsets of pa 5- 15" (w 30-75 kpc). In the former scenario, 
unless the galaxies have L ~ 0.1 L*, the resulting redshift dis- 
tribution is bimodal with peaks at z ~ 0.5 and ~ 3. Such a 
scenario most likely requires a bimodal age distribution in the 
context of NS-NS mergers, or two distinct progenitor systems 
dominating at low and high redshifts. While this cannot be 

8 We ignore th e fac tor of ~ 5 correction to the energy, due to incomplete 
energy coverage ( 35. 2t . since this would roughly apply to all bursts in roughly 
the same way. The higher redshift events may require somewhat smaller 
corrections due to redshifting of the spectral peak into the softer energy band 
covered by Swift. 



ruled out by present data, the lack of overlap at z ~ 1 -2 is dif- 
ficult to explain. We also conclude that it is unlikely that any 
faint underlying hosts are low-luminosity galaxies (L < 0.01 
L*) at similar redshifts to the detected hosts, since this does 
not naturally explain the difference in afterglow and prompt 
emission properties. 

In the context of the large offsets scenario, the probability 
of chance coincidence (both individually and for the sample) 
is lower for the brighter galaxies at offsets of ss 5- 15" than 
for the faint galaxies at offsets of « 2". However, it is not 
unlikely that the true associations are a mix of both popula- 
tions, since in each case there is a non-negligible probability 
that 1-2 of the associations are spurious. We note that optical 
positions are crucial for a uniform comparison of the chance 
coincidence probabilities for faint galaxies with small offsets 
and bright galaxies with larger offsets; in the case of only X- 
ray positions there is an inherent bias (in the sense of a poste- 
riori) against associations with faint galaxies inside the error 
circle, despite the fact that they are consistent with no offset. 

From spectroscopic observations for 3 of the 5 bursts 
we find that the galaxies with the lowest probability of 
chance coincidence are a star forming galaxy at z = 
0.111 (GRB 061201), an early-type galaxy at z = 0.473 
(GRB 070809), and an early-type cluster member galaxy at 
z = 0.403 (GRB 090515). If these associations are indeed cor- 
rect, they only slightly alter the host demographics, which are 
dominated by star forming galaxies (Berger 2009). However, 
this does suggest that our present understanding of the relative 
ratios of star forming and elliptical hosts may be incomplete. 

The resulting distributions of angular, physical, and host- 
normalized offsets for the bursts with and without coincident 
hosts appear to be distinct, rather than continuous. How- 
ever, the joint distribution of projected physical offsets is in 
good agreement with theoretical predictions for NS-NS bi- 
nary mergers. On the other hand, the predicted distribution for 
dynamically-formed NS-NS binaries in globular clusters pro- 
vides a much poorer fit to the entire data set, unless they ac- 
count for only the bursts with large offsets (Sample 2). In the 
case of a hybrid population of primordial and dynamically- 
formed binaries, with the latter accounting for only the large 
offsets, the fact that all the remaining offsets (Sample 7) are 
< 10 kpc, is indicative of no significant kicks. The large phys- 
ical offsets also naturally explain the fainter afterglow emis- 
sion as a result of lower circumburst densities. The resulting 
isotropic 7-ray energies match the observed distribution for 
short GRBs with coincident hosts, either at z ~ 0.3 or at z > 1 . 

Our conclusion that large offsets of ~ 15-70 kpc (corre- 
sponding to rs 4— 10 galactic effective radii) are a likely ex- 
planation for the bursts with optical afterglows and no coin- 



cident hosts is of fundamental importance. The only progen- 
itor model that naturally explains this result in the merger of 
NS-NS/NS-BH binaries, most likely due to kicks, or possibly 
with a minor contribution from a globular cluster population 
(accounting specifically for the events in Sample 2). While 
a conclusive demonstration of a large offset requires an ab- 
sorption redshift measurement that matches an offset galaxy 
redshift, the distribution of optical afterglow magnitudes in- 
dicates that this will be difficult to achieve. Indeed, an ab- 
sorption redshift is a vailable for only one likely short GRB 
JLevesqueetal.ll2010h . 

We end with several important observations. First, rapid op- 
tical follow-up of short GRBs with 8-meter class telescopes is 
essential since observations to a depth of about 25 mag within 
a few hours after the burst may recover nearly all optical after- 
glows, regardless of circumburst density (Figure |5). Second, 
short GRBs with only X-ray positions are unlikely to provide 
strong support for either negligible or large offsets due to the 
appreciable size of the error circles and the fact that the sam- 
ple with optical afterglows exhibits wide host galaxy diversity, 
i.e., events with coincident bright hosts at z ~ 0.2- 1, events 
with coincident faint hosts likely at z > 1, and events with 
no coincident hosts likely due to offsets. Third, a meaningful 
study of short GRB offsets requires a statistical approach to 
mitigate the shortcomings of a posteriori chance coincidence 
probabilities, as well as to incorporate the relevant informa- 
tion from afterglow and prompt emission observations. 

Finally, we note that from a wide range of observations of 
both the afterglows and host galaxies it appears that the case 
for NS-NS mergers as the progenitors of short GRBs is gain- 
ing observational support. Our main result here is that short 
GRBs with optical afterglows and no detected host galaxies 
are somewhat less likely to be explained with high redshifts or 
with dwarf galaxy hosts at low redshift. Instead, these bursts 
likely exploded ~ 15 kpc from galaxies at z ~ 1 or tens of kpc 
from galaxies at z ~ 0.3. With larger samples, these possibil- 
ities will allow us to improve our understanding of the short 
GRB redshift distribution, the host galaxy demographics, and 
predictions for gravitational wave detections. If high redshifts 
indeed turn out to be prevalent, this will have a significant 
effect on the possibility of multiple progenitor populations. 
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TABLE 1 
Properties of Short GRBs with Optical Afterglows or Limits 
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GRB 


I'm) 


z a 


F h 


tx 


Fx c 


fopt 


^V.opt 


Refs. 




(s) 




(erg crrT 2 ) 


(hr) 


(erg crrT- s ) 


(hr) 


(A«Jy) 






Short GRBs with 


Optical Detections 


(Samples 1 & 2) 






050709 


0.07 


0.161 


2.9 x 10-' 


60.5 


3.5 x 10" 15 


34 


2.3 


1-3 


050724 


3.0 


0.257 


3.9 x 10~ 7 


11.8 


1.7 x 10~ 12 


12 


8.4 


4-6 


051221A 


1.4 


0.546 


1.2 x 10- 6 


3.1 


1.8 x U)- n 


3.1 


5.8 


7-8 


060121 


2.0 


<4.3 


4.7 x 10- 6 


6.1 


1.8 x 10-' 2 


7.4 


8.8 


9-11 


060313 


0.70 


<3 


1.1 x 10~ 6 


2.8 


4.0 x 10~ 12 


2.8 


10.8 


11-12 


061006 


0.42 


0.438 


1.4 x 10- 6 


19.5 


1.0 x U)- 13 


14.9 


2.9 


11,13 


061201°' 


0.80 


< 1.7, 0.111' 


? 3.3 x 10~ 7 


8.6 


2.5 x 10- 13 


8.6 


2.9 


14-15 


070707 


1.1 


<4.3 


1.4 x 10~ 6 


12.5 


2.3 x 10~ 13 


11 


1.9 


16 


070714B 


3.0 


0.923 


7.2 x 10~ 7 


24.4 


1.5 x 10~ 14 


23.6 


0.7 


17 


070724 


0.40 


0.457 


3.0 x 10~ 8 


2.35 


5.5 x 10- 13 


2.3 


5.0 


18-19 


070809 d 


1.3 


< 3, 0.473? 


1.0 x 10~ 7 


11.0 


3.0 x 10~ 13 


11 


0.8 


20-21 


071227 


1.8 


0.381 


2.2 x 10~ 7 


7.0 


6.9 x 10~ 14 


7.0 


1.6 


13 


080503 d 


0.32 


<3, ■■■ 


6.1 x 10~ 8 


66.0 


< 7.8 x 10~ 14 


25.9 


0.3 


22 


080905 


1.0 


0.122 


1.4 x 10~ 7 


18.4 


< 5.9 x 10- 14 


8.5 


0.8 


23 


090305 d 


0.40 


<4.3, ■■■ 


7.5 x 10~ 8 


0.45 


< 8.5 x 10~ 14 


0.45 


1.5 


24-25 


090426 


1.28 


2.609 


2.5 x 10~ 7 


2.0 


1.1 x 10~ 12 


2.6 


20.0 


26-27 


090510 


0.30 


0.903 


3.4 x 10~ 7 


6.4 


1.7 x icr 13 


9.0 


2.3 


28 


09051 5 d 


0.04 


< 4.3, 0.403' 


.' 2.1 x 10~ 8 


4.9 


< 9.2 x 10~ 14 


1.9 


0.1 


29 


091109B 


0.30 




1.9 x 10~ 7 


3.8 


1.9 x 10~ 13 


5.7 


0.5 


30 


100117 


0.30 


0.92 


9.3 x 10~ 8 


8.4 


< 2.0 x 10-' 4 


8.4 


0.3 


31 




Short GRBs 


with Optical Limits (Sample 3) 








050509B 


0.04 


0.225 


9.5 x 10~* 


9.1 


< 2.0 x 10~ 14 


2.1 


0.7 


32-33 


050813 


0.60 




1.2 x 10~ 7 


94.9 


< 2.6 x 10-' 4 


12.8 


1.9 


34-36 


051210 


1.27 




8.1 x 10~ 8 


6.9 


< 2.9 x 10~ 14 


19.2 


1.6 


11,37 


060502B 


0.09 




4.0 x 10~ 8 


18.6 


<9.1 x 10~ 14 


16.8 


0.7 


11,38 


060801 


0.50 


1.130 


8.1 x 10~ 8 


10.9 


< 9.0 x 10~ 15 


12.4 


0.8 


11 


061210 


0.19 


0.409 


1.1 x 10~ 6 


72.4 


< 6.7 x 10~ 14 


2.1 


1.4 


11 


061217 


0.21 


0.827 


4.6 x 10~ 8 


38.6 


< 1.7 x 10~ 14 


2.8 


2.0 


11 


070429B 


0.50 


0.902 


6.3 x 10~ 8 


11.8 


< 3.8 x 10~ 14 


4.8 


0.6 


39 


080426 


1.30 




3.7 x 10~ 7 


6.1 


2.9 x 1(T 13 


7.5 


2.6 


40 


080702 


0.50 




3.6 x 10~ 8 


2.9 


< 4.0 x 10~ 14 


12.1 


12. 


41 


081226 


0.40 




9.9 x 10~ 8 


3.2 


< 2.9 x 10-' 4 


1.6 


0.9 


42 


100206 


0.12 




1.4 x 10~ 7 


6.9 


< 2.0 x 10~ 14 


15.7 


0.5 


43 



NOTE. — Prompt emission and afterglow data for short GRBs with detected optical afterglows (top section) and 
deep optical afterglow limits (bottom section). 

a Redshifts include spectroscopic measurements, limits from afterglow detections in the UV/optical, and for the 
bursts in Sample 2, redshifts for galaxies with the lowest probability of chance coincidence (marked by ?). 
'' The fluences are in the observed 15-150 keV band, with the exception of GRB 050709 (2-400 keV) and 

GRB 060 1 2 1 (2 - 400 keV). 

c All XRT data are from lEvans et all 120071) and lEvans etal] 120091) . 

Short GRBs in Sam ple 2. 

Refere nces: [11 Villasenor etai] I2005T); [21 IFox et all 120051); 131 IHiorth eTaT] I200 5[): T41 IBarthelmv et all 
12005); [5] Bergeretal. (2005); [6] Grupe et al. i 2006); [7] Burrows et al. (2006); [8] Soderb erg et all 120061) ; 
[9] de Ugaite Postigo et al. (2006); [10] Levan et al. 1 2006a); [11] Bergeret al. (2007); [12] Roming et al. (2006); 
[131 ID'Avanzo et all 120091) ; 114 1 IStratta et all 120071) ; TlSl IFong et all 120101); T16 1 Piranomo nteetal] |2008[) ; 
[171 IGraham et all (20 09); [18] Berger et al \ j2 009t); 1191 IKocevski et all 1 2010); [20] Perlev et all 1200% ; 
r211 iPerlev et all (20081) ; f221 iPerlev etal . i2009); [23] Row linson et all I2O10I); [241 ICenko et a l. i 2009); [25] 
Berger & Kelson (2009); [26] Antonelli et al. (2009); [27] Levesque et al. .2010); [28] McBr een et al] 12010) ; 
[29] Rowlinson et al. (2010); [30] Malesani et al. (2009); [31] Fon g et al in pre p.; [321 [ Gehrels et alj 120051) ; 
[331 JBloom eta l. (2006); [34] Fe rrero et al. i 2007); [35] Berger 1 2006b); [36] Proch aska et alj 120061); [37 ] 
La Parola et al. (2006); [38] Bloom et al. (2007); [39] Cenko et al. (2008); [401 Ide Ugarte Postigo et all 120081) ; 
[411 IGrecoet alj (20081) ; [421 1Berger etal. 1 2008); [43] Berger et al|(201CI) 
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TABLE 2 

Observations of Short GRBs with Optical 

Afterglows and no Coincident Host Galaxies 

(Sample 2) 



GRB 


Instrument 


Filter 


(s) 


(AB mag) 


061201 


HST/ACS 


F814W 


2224 


26.0 


070809 


Magellan/LDSS3 


r 


1500 


25.4 


080503 


HST/WFPC2 


F606W 


4000 


25.7 


090305 


Magellan/LDSS3 


r 


2400 


25.6 


090515 


Gemini-N/GMOS 


r 


1800 


26.5 



Note. — a Limits are 3<t. 
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FIG. 1. — ffi>T/ACS/F814W image of the location of GRB 061201. Galaxies near the position of the optical afterglow (cross-hairs) are marked. 
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FIG. 2. — Magellan/LDSS3 r-band images of the location of GRB 070809. Galaxies near the position of the optical afterglow (cross-hairs) are marked. 
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Fig. 3. — //577WFPC2/F606W image of the location of GRB 080503. Galaxies near the position of the optical afterglow (cross-hairs) are marked. A faint 
galaxy at a separation of only 0.8" was found by Perley et al. (2009) based on a deeper stack of HST/WFPC2 observations. These authors also find that the 
galaxy marked "S5" is located at z = 0.561, leading to a physical offset of 85 kpc. 
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FIG. 4. — Magellan/LDSS3 r-band images of the location of GRB 090305. Galaxies near the position of the optical afterglow (cross-hairs) are marked. 
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FIG. 5. — Gemini-North/GMOS r-band images of the location of GRB 090515. Galaxies near the position of the optical afterglow (cross-hairs) are marked. 
Note that the object coincident with the cross-hairs is the optical afterglow. 
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FIG. 6. — Magellan/LDSS3 spectrum of the galaxy with the lowest probability of chan ce coincidence near the position of GRB 061201. T his galaxy is marked 
"S4" in Figure[T] It has aredshift of z = 0.1 11 and it is undergoing active star formation iBerger 2006a; Stratta et al. 2007; Fong et al. 2010). 
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FIG. 7. — Magellan/LDSS3 spectrum of the galaxy with the lowest probability of chance coincidence near the position of GRB 070809. This galaxy is marked 
"S3" in Figure[2] It has a redshift of z = 0.473 and is an early-type galaxy with no evidence for on-going star formation activity (see inset). 
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FIG. 8. — Magellan/LDSS3 spectra of three galaxies with a low probability of chance coincidence near the position of GRB 090515. The galaxy with the lowest 
probability of chance coincidence is marked "S5" in Figure|5] It has a redshift of z = 0.403 and is an early-type galaxy which is part of a galaxy cluster (Fong et 
al. in prep.) 
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FIG. 9. — Optical afterglow brightness on timescales of a few hours after the burst for short GRBs with detected afterglows (Sample 1: black squares; Sample 
2: red squares) or upper limits (gray triangles). The lines at the top right indicate the fading tracks for afterglow decay rates of a = —0.5 and — 1. The right panel 
shows the projected histogram for the bursts with detected afterglows (hatched) and upper limits (open). The symbols mark the mean for each sample, and the 
vertical bar marks the standard deviation for Sample 1 . 
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FIG. 10. — Cumulative afterglow brightness distributions for the bursts in Figure[9] extrapolated to a common fiducial time of 8 hr after the burst with a fading 
rate of a = —0.75. The K-S probabilities relative to the sample with detected afterglows and coincident hosts are noted in the figure. It appears unlikely that the 
bursts with no coincident hosts, and those with deep upper limits, are drawn from the same distribution as the bursts with detected hosts. 
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FIG. 1 1 . — Histograms of 7-ray fluence (top), afterglow X-ray flux at 8 hr (middle), and duration (bottom) for the three short GRB samples discussed in this 
paper. The arrows mark the mean for each sample, indicating that the bursts in Sample 2 and Sample 3 have lower 7-ray fluences, fainter X-ray fluxes, and shorter 
durations, than the bursts with detected afterglows and coincident hosts. 
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FIG. 12. — Short GRB 7-ray fluence as a function of duration for the three samples discussed in this paper. An overall correlation is apparent in the data. The 
bursts in Sample 2 and Sample 3 appear to lie below the mean correlation for the bursts in Sample 1, i.e., they have lower fluences for their durations, or longer 
durations for their fluences. 
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FIG. 13. — X-ray versus optical flux for the bursts from Figure[9] The cross-hatched region marks the median and standard deviation assuming the expected 
power-law correlation with an index /3ox- The light shaded region marks the region occupied by long GRBs I Jakobsson et al. 2004). The distributions for long 
and short GRBs are largely indistinguishable, as are the distributions for short GRBs with an without coincident hosts. We note that a large fraction of the bursts 
with optical upper limits also have undetected X-ray afterglows on timescales of a few hours after the burst. The overall similarity between the ratio of optical to 
X-ray flux for long and short GRBs does not allow us to clearly locate the synchrotron cooling frequency (u c ) in relation to the X-ray band. If v c > mix for short 
GRBs, the faintness of the optical afterglows for bursts with no coincident hosts cannot be used to distinguish density and redshift effects. 
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FIG. 14. — Probability of chance coinc idenc e as a function of distance from a short GRB optical afterglow position for galaxies near the location of each burst. 
These are the galaxies marked in Figures [TT5l In each panel we mark the galaxy with the lowest probability of chance detection with a circle. In 4 of the 5 cases, 
the lowest probability is associated with galaxies that are offset by ~ 5- 15". Moreover, even the nearest galaxies are offset by PS 1.6-5.8". 
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FIG . 15. — Same as Figure [141 but for short GRBs with coincident faint hosts. In this case, the lowest probability of chance coincidence is associated with the 
underlying faint host. 
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FIG. 16. — Host galaxy optical magnitude as a function of redshift for short GRB hosts (black squares), long GRB hosts (gray circles), and galaxies with a 
luminosity of 0. 1 - 1 L* (shaded region). The dashed lines mark the upper limits at the GRB positions for the short GRBs with no coincident hosts. The arrows 
mark the upper limits on the redshifts of three bursts with faint hosts, based on the detection of the afterglows in the optical band (i.e., lack of a strong Lyman 
break). If underlying host galaxies exist for Sample 2, their non-detection indicates z > 1 .5 (for 0. 1 L* ) or > 3 (for L* ). The alternative possibility that they are 
located at similar redshifts to the detected hosts, requires < 0.01 L*, but this does not naturally explain their fainter afterglows. 
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FIG. 17. — Histogram of projected angular offsets relative to the host galaxy center for short GRBs with coincident hosts (hatched black), and bursts with no 
coincident hosts if the galax ies w ith lowest chance coincidence probability are the hosts (hatched red), or if the faint galaxies with smallest angular separation are 
hosts (open red); see Figure [\A\ The dashed line is a log-normal fit to the bursts with coincident hosts. 
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FIG. 18. — Same as Figure fTTI but normalized relative to the host effective radii, R e . The dashed line is a log-normal fit to the bursts with coincident hosts. 
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FIG. 19. — Same as Figure [TTl but projected physical offsets in units of kpc. The dashed line is a log-normal fit to the bursts with coincident hosts. 
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FIG. 20. — Cumulative distributions of projected physical offsets for short GRBs with coincident hosts (black line), and combined with offsets for the hosts 
with the lowest probabilities of chance coincidence (thick red line) or the faint hosts with smallest angular offsets (thin red line). Also shown are predicted 
distributions for NS-NS kicks from several models (Bloom et al. 1999; Fryer et al. 1999; Belczynski et al. 2006), and for dynamically-formed NS-NS binaries 
from globular clusters (shaded region marks a range of predictions for host galaxy masses of 5 X 10 10 — 10 12 Mq; Salvaterra et al. 2010). The models with kick 
velocities are in good agreement with the measured offset distribution for either set of galaxy associations, while the globular clusters model provides a poor 
match to the data. 
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FIG. 21. — Isotropic-equivalent 7-ray energy as a function of redshift for the bursts with detected optical afterglows. We plot the inferred energies for Sample 
2 at the redshifts corresponding to the galaxies with the lowest probability of chance coincidence (red squares); at z ~ 1 corresponding to the faint galaxies with 
the smallest angular separations (thick open squares); and at z ~ 3 corresponding to the case of undetected underlying hosts with a luminosity of about L* (thin 
open squares). The bursts in Sample 2 match the distribution for the Sample 1 bursts at any of these redshift intervals. 



